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A.  Summary 

To  date,  essentially  all  thermal  and  thermoelectric  materials  and  devices  are  based  on  linear- 
response  phenomena,  where  the  effect  (for  example  flow  of  heat  or  charge)  is  proportional  to  the 
driving  force  (a  thermal  or  electric  bias).  This  is  in  contrast  to  electronic  devices,  where  the  use 
of  nonlinear  elements,  such  as  diodes,  is  widespread  and  enables  key  functionalities. 

This  report  addresses  thermal  phenomena  in  nonlinear  response,  with  the  goal  to  explore 
fundamentally  new  effects  and  functionalities  that  may  enable  game-changing,  novel  devices, 
and  enhanced  performance. 
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Electrons  in  semiconductor  devices  with  sufficiently  small  (submicron)  features  move  on  straight 
trajectories  (ballistically)  between  device  boundaries,  and  electronic  properties  are  therefore 
determined  by  the  device  symmetry.  In  this  highly  controllable  system  we  have  observed 
transverse  thermoelectric  effects  in  multi-terminal  devices,  and  we  have  developed  a  scattering- 
matrix  theory  that  explains  these  effects.  Specifically,  using  this  theory,  we  can  predict  the 
magnitude  and  direction  of  the  transverse  thermoelectric  effect  based  on  (nonlinear)  electric 
characterization  of  the  device.  In  other  words,  we  now  understand  this  system  sufficiently  well 
that  we  can  predict  its  response  to  a  thermal  bias  from  measurements  using  an  electric  bias  alone. 
We  have  thus  established  a  framework  for  the  targeted  development  of  novel  functional  materials 
and  devices  based  on  nonlinear  thermal  phenomena. 

The  insights  gained  may  be  applicable  not  only  to  electronic  systems,  but  also  to  heat  flow 
carried  by  phonons  in  thermal  devices.  Military  applications  of  rectified  heat  flow  include  solid- 
state  refrigerators  for  FIR  detectors  (reduction  of  parasitic  back  flow  of  heat),  and  thermal 
management  in  high-power  electronics. 


B.  Introduction 

B.l.  Background:  New  functionality  in  the  nonlinear  response  regime 

The  functionalities  of  essentially  all  current  applications  of  thermoelectricity  and  of  thermal 
management  are  well  described  within  linear  response  theory:  Ohm’s  Law  (7  =  GV)  for  charge 
current,  Fourier’s  Law  ( dQ/dt  =  KAT)  for  thermal  conduction. 

In  the  nonlinear  regime,  the  higher  terms  in  the  following  expansions  become  important: 

I=GV+G2V2  +  G3V3  +  ... 
dQ/dt  =  KAT+K2Af  +  K3Af 

Current  significant  progress  in  nanostructured  materials  offers  strong  motivation  to  explore 
nonlinear  thermal  phenomena.  First,  on  length  scales  shorter  than  the  characteristic  energy 
relaxation  lengths  of  electrons  and/or  phonons,  nonlinear  effects  (the  breakdown  of  Fourier’s 
Law)  are  very  likely  because  the  energy  of  electrons  and/or  phonons  is  no  longer  controlled  by 
the  temperature  of  the  bath,  but  by  the  applied  bias.  Second,  the  ability  to  structure  materials  on 
the  nanoscale  with  very  high  precision  can  be  used  to  control  the  material’s  symmetry,  enabling 
the  control  of,  for  example,  rectifying  phenomena. 

New  functionality  is  possible  in  this  regime.  For  example,  the  second  term  in  each  equation  above 
signifies  rectifying  effects,  because  this  term  enhances  the  current  for  one  sign  of  the  bias,  and 
reduces  it  for  the  opposite  sign. 

B.2.  Motivation:  Enabling  novel  functionalities  and  enhanced  performance 

The  ability  to  control  nonlinear  thermal  behavior  would  enable  a  number  of  novel  functionalities 
and  applications,  as  illustrated  in  the  following  examples: 

Example  1:  Rectified  heat  flow  in  Peltier  coolers,  or  in  thermal  management.  In  Peltier 
coolers,  an  applied  electric  current  removes  heat  from  a  material,  such  that  it  can  be 
cooled  below  ambient  temperature.  The  performance  of  Peltier  coolers  is  limited  by  the 
backflow  of  heat,  in  a  direction  opposite  to  that  of  the  electric  current.  A  nonlinear 
material  in  which  electrons  move  easily  in  the  direction  of  the  current,  but  have  a  low 
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mobility  in  the  direction  of  the  heat  backflow,  could  have  significantly  enhanced 
performance.  To  date,  thermal  rectification  has  been  demonstrated  at  some 
heterostructure  interfaces  but,  to  our  knowledge,  no  functional,  nanostructured  material 
with  this  property  has  been  designed. 

Example  2:  Transverse  thermoelectrics.  Also  the  performance  of  thermoelectric  power 
generators  is  limited  by  undesired  heat  flow,  but  in  this  case  electric  current  and  parasitic 
heat  flow  are  in  the  same  directions.  One  relatively  unexplored  possibility  is  to 
geometrically  decouple  charge  and  heat  flow,  by  using  transverse  thermoelectric  effects, 
where  the  charge  flow  is  perpendicular  to  the  applied  thermal  gradient.  Such  effects  can 
be  induced  in  asymmetric  nanostructures,  potentially  enabling  the  design  of  functional  TE 
materials  with  enhanced  performance. 

Here,  we  report  our  research  progress  on  the  physical  mechanisms  underlying  nonlinear  thermal 
effects  in  a  highly  controlled  experimental  system  that  is  amenable  to  the  development  and  test  of 
a  rigorous  theoretical  framework. 


C.  Methods,  Assumptions,  and  Procedures 

C.l.  Device  Fabrication 

Our  experiments  were  carried  out  on  two-dimensional  electron  gas  (2DEG)  devices  embedded  in 
InGaAs/InP  heterostructure  wafers.  These  heterostructures  are  grown  on  the  surface  of  InP  wafers 
at  Lund  University  using  metal  organic  vapor  phase  epitaxy.  The  composition  of  the 
heterostructure  is  shown  in  Fig.  1(a),  the  details  of  which  are  covered  in  Refs.  2-4.  Wafers  grown 
in  this  material  and  manner  are  nearly  perfectly  flat,  offer  high  electron  mobilities  (pe  >  105 
cm2/Vs) 2,3 ,  and  can  be  chemically  etched  without  the  risk  of  surface  oxidization  when  exposed  to 
air.  All  fabrication  steps  used  to  create  the  device  geometry  on  the  2DEG  wafer  are  done  using 
electron  beam  lithography  (EBL),  wet  chemical  etching,  metallization,  and  rapid  thermal 
annealing.  The  fabrication  details  can  be  found  in  Ref.  5. 


C.2.  Measurement  Methods 

Two  scanning  electron  micrographs  of  the  ballistic  electron  devices  used  in  our  experiments  are 
shown  in  Fig.  l(b-c),  labeled  as  device  1  and  2  respectively  from  here  on.  Both  figures  show  the 
intersection  of  four  leads  (terminals),  containing  an  asymmetric  (triangular)  antidot,  a  scatterer 
for  electrons.  An  important  property  of  device  2  is  the  complete  symmetry  of  the  device,  except 
the  central  triangular  scatterer.  The  electrons  in  the  terminals  are  heated  using  two  opposing,  180° 
out  of  phase,  37  Hz  AC  currents,  v±H,  applied  to  the  individual  heating  channels.  Fig.  l(b-c)  show 
example  circuit  diagrams  for  these  measurements.  The  lco  (37  Hz)  voltage  at  the  heated  terminal, 
as  measured  by  the  voltage  probes  Pi  and  Vp, 2,3,4)  in  the  respective  devices,  is  minimized  by 
tuning  the  relative  amplitudes  of  v±H.  Any  DC  voltage  measured  here  is  also  removed  by 
applying  a  DC  shift,  Vue,  to  v±H.  The  thermovoltage  responses  are  then  measured  by  measuring 
the  voltages  at  each  terminal  using  the  available  voltage  probes.  Using  this  method,  we  can 
measure  the  thermovoltage  response  due  to  heating  the  electrons  in  terminals  1  or  3  for  the  device 
1,  and  any  of  the  four  terminals  for  the  device  2.  When  heating  terminal  3  in  device  1,  however, 
we  cannot  minimize  the  DC  and  loo  voltage  due  to  the  lack  of  a  side  voltage  probe. 
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The  junctions’  response  to  an  electronic  excitation  was  also  measured  to  characterize  their  linear 
and  nonlinear  transport  properties.  The  electronic  measurements  were  conducted  by  applying  a 
controlled  electric  current  between  two  terminals  while  measuring  the  voltage  responses  at  each 
terminal.  Due  to  the  geometry  of  device  1,  the  device  could  only  be  electrically  biased  between 
terminals  1  and  3.  See  Refs.  1,  5,  and  7  for  additional  details  on  both  measurements. 

In  addition  to  applying  the  electrical  and  thermal  biases  to  device  2,  we  could  also  apply  a 
magnetic  field  perpendicular  to  the  plane  of  the  2DEG  to  study  the  quantum  thermoelectric 
transport  symmetries.  In  addition  to  the  varying  electrical  and  thermal  excitation  stated  above,  our 
symmetry  measurements  consist  of  a  constant  electrical  or  thermal  excitation  with  a  varying 
magnetic  field.  For  each  of  these  symmetry  measurements,  multiple  four-terminal  combinations 
are  used  to  gather  the  required  information  for  analyzing  the  junction’s  symmetry. 


Fig.  1.  In  (a),  a  cross-section  of  the  heterostructure  InGaAs/InP  wafer  is  shown;  the  dopant  layer  is  depicted 
by  the  black  band  and  the  2DEG  the  green  band.  In  (b),  a  typical  measurement  setup  for  a  four-terminal 
junction  is  shown  around  a  scanning  electron  micrograph  (SEM)  of  a  junction.  In  (c),  a  second  device  is 
shown  that  has  been  tailored  to  study  the  magnetic-field  symmetries  of  quantum  thermoelectric  transport  in 
these  four-terminal  junctions.  The  central  junction  in  (c)  is  similar  in  structure  and  size  to  the  junction 
shown  in  (b).  In  both  (b)  and  (c)  the  four  terminals  in  the  junction  are  numbered  1  through  4.  The  white 
circles  on  the  edges  of  the  two  devices  indicate  the  contact  points  to  the  side  channels,  which  are  used  to 
bias  the  devices  both  electrically  and  thermally.  Panel  (b)  was  taken  from  Ref.  1 ,  and  panel  (c)  was  taken 

from  Ref.  7. 
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D.  Results  and  Discussion 

D.l.  Transverse  Thermovoltage  Response 

The  key  experiment  is  to  apply  a  heating  current  to  terminal  1 ,  and  to  measure  the  transverse 
voltage  response  across  terminals  2  and  4.  The  first  four  harmonics,  with  the  first  being  the 
frequency  of  the  heating  voltage,  of  the  transverse  thermovoltage  show  that  the  2nd  harmonic 
dominants  the  transverse  response,  Fig.  2(a).  Since  Joule  heating  generates  the  temperature 
increase  in  the  terminal,  a  dominant  2nd  harmonic  is  indicative  of  a  linear  response  to  the  applied 
temperature  differential. 

The  literature6  also  shows  a  dominant  2nd  harmonic  due  to  electrical  excitation.  We  were  able  to 
confirm  that  the  transverse  response  we  observe  is  due  to  thermal  effects  by  measuring  the  lateral 
voltage  drop  across  the  junction  during  the  heating  measurements,  Fig.  2(a)  inset,  and  comparing 
it  to  the  junction’s  response  to  electrical  excitations,  Fig.  2(b).  We  find  that  the  resulting  2co 
response  due  to  the  lateral  voltage  drop  is  around  10  nV,  which  is  roughly  100  times  smaller  than 
the  observed  thermoelectric  response. 


Fig.  2.  Panel  (a)  shows  the  first  four  harmonics  of  the  transverse  thermovoltage.  The  dominance  of  the  2co  signal 
indicates  that  the  response  in  linear  in  the  applied  temperature  differential.  The  inset  shows  the  measured  lateral 
voltage  drop  during  the  applied  thermal  gradient  measurements.  In  panel  (b)  we  checked  that  the  2co  observed 
response  is  not  due  to  electronic  excitations  by  characterizing  the  junction’s  electronic  response.  Figures  taken  from 

Ref.  1. 


D.2.  Observed  Left-Right  Asymmetry  and  the  Two  Quantum  Point  Contact  Model 

Based  on  the  left-right  geometric  symmetry  of  the  junction,  a  naive  guess  would  say  that  the 
transverse  thermovoltage  should  not  change  upon  reversing  the  lateral  thermal  gradient,  that  is, 
V24(0i,03)  «  V24(03,0i).  Upon  reversing  the  lateral  thermal  gradient  however,  we  find  that 
transverse  thermovoltage  changes  sign. 

This  observation  can  be  understood  by  conceptually  treating  the  junction  as  two  quantum  point 
contact  (QPC)  pairs.  The  thermopowers  of  the  top-right  and  bottom-right  QPCs  determine  the 
junction’s  response  to  heating  terminal  1.  Similarly,  the  thermopowers  of  the  left  QPC  pair 
determines  the  response  to  heating  terminal  3.  Since  a  QPC’s  thermopower  is  a  series  of  peaks  as 
a  function  of  the  Fermi  energy  (occupation  umber),  the  sign  change  in  our  device  can  easily  arise 
due  to  a  slight  left-right  asymmetry  in  the  QPC  widths. 

The  conceptual  two  QPC  model  is  also  able  to  give  a  simplistic  explanation  about  the  non- 
linearities  in  the  thermovoltages,  indicated  by  the  vertical  dashed  lines  in  Fig.  3.  As  the 
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temperature  drop  across  a  QPC  increases,  its  total  thermopower  increases  due  to  thermal 
broadening  about  the  Fermi  energy  accessing  neighboring  thermopower  peaks.  The  downward 
change  in  slope  in  the  heat-right  curve  can  then  be  interpreted  as  an  increase  in  the  number  of 
accessed  thermopower  peaks  in  the  top  right  QPC  as  a  function  of  0i.  Similarly  for  the  heat-left 
configuration,  the  two  nonlinear  points  indicate  increases  in  the  number  of  accessed  thermopower 
peaks  in  the  bottom  left  QPC. 


Fig.  3.  Upon  reversing  the  lateral  thermal  gradient,  we  see  that  the  transverse  thermovoltage  changes  sign, 
something  initially  unexpected  based  on  the  geometrical  symmetry  of  the  junction.  The  vertical  dashed  lines 
roughly  indicate  the  onset  of  the  non-linear  regions  in  the  transverse  responses.  Figure  taken  from  Ref.  1. 


D.3.  Multi-Terminal  Thermoelectric  Scattering  Theory 

Though  the  two  QPC  model  is  a  useful  model,  all  four  QPCs  in  the  junction  interact  with  each 
other,  something  the  two  QPC  model  neglects.  To  fully  treat  the  four-terminal  junction,  we 
developed  a  multi-terminal  scattering  theory  capable  of  describing  any  ballistic,  mesoscopic 
junction  up  to  quadratic  order  in  the  voltages  and  linear  order  in  the  temperatures.  We  found  that 
the  current  through  the  ath  junction,  Ia,  can  be  described  as 


P 


G«p(Va  ~  Vf)  +  Efc  -  Vr)Gc.l>r(Va  ~  Vp)+L«v(0"  ~  9P) 


2e 

Ga(j  =  tap  C &F  ) 


G  e3  dap 


aPr 


h  dE 


2e2  3, 


}Pr' 


ap 


h  dV„ 


L 


•ap 


iTpekjO  dtap 


3  h  dE 


where  e  is  the  electron  charge,  h  is  Planck’s  constant,  kp  is  the  Boltzmann  constant,  0 is  the 
background  temperature,  Ep  is  the  common  Fermi  energy,  ta/3  is  the  transmission  probability 

between  the  ath  and  pth  terminals,  Va  is  the  voltage  in  the  ath  terminal,  6a\ s  the  voltage  in  the  ath 
terminal. 
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Using  this  scattering  theory,  we  were  able  to  estimate  the  temperature  increase  in  terminal  1  as  a 
function  of  the  applied  heating  voltage  and  thus  estimate  the  effective  transverse  thermopower  of 
the  junction1.  We  found  a  transverse  thermopower  of  roughly  600  nV/K,  which  is  comparable  to 
the  difference  between  two  QPC  thermopowers. 

D.4.  Magnetic-Field  Symmetries  in  Quantum  Thermoelectric  Transport 


Based  on  micro-reversibility,  the  thermoelectric 
transport  properties  of  our  devices  should  be 
symmetric  under  a  reversal  of  the  magnetic  field 
if  the  current  and  voltage  leads  are  exchanged  at 
the  same  time.  For  electrically  biased  transport, 
this  leads  to  specific  symmetries  in  the 
transmission  function,  tap  ( E,B )  =  t pa  ( E-B ), 
where  tap  is  the  transmission  function  between 

terminals  a  and  p,  E  is  the  electron  energy,  and  B 
is  the  magnetic  field.  These  symmetries  have  been 
thoroughly  investigated  in  past  literature  through 
measurements  using  electronic  excitation.  Since 
the  thermoelectric  coefficients,  Lap,  are 
dependent  on  the  energy  derivative  of  the 
transmission  function,  dEtap,  we  expect  that  the 
same  symmetries  should  be  present  in  the 
thermally  driven  case, 

Laf(B)  =  Lfa(rB). 

Our  observations  show  that  both  the  electronic 
and  thermal  linear-symmetries  are  present  in  our 
devices,  Fig.  4.  Figure  4(a)  shows  two  expected 
symmetries  based  on  the  measured  four-terminal 
electrical  resistances,  Rappy  =  Vap  /  hy  the 

symmetries  of  which  are  RapSy{B )  =  RSrap  (-5) . 

The  newly  observed  thermal  bias  symmetries  are 
shown  in  Fig.  4(b-d). 

Though  the  magnetic  field  traces  in  Fig.  4  show 
the  expected  symmetries,  the  symmetries  are  not 
as  perfect  as  expected.  Significant  small-scale 
fluctuations  are  present  in  the  electric  and  thermal 
bias  measurements,  though  smaller  in  relative 
magnitude  for  the  former.  Two  possible 
explanations  for  the  deviation  in  the  thermal  bias 
measurements  are:  i)  unintentional  heating  of  the 
three  non-heated  terminals,  for  example  terminals 
2,  3,  and  4  in  Fig.  1(c);  and  ii)  inelastic  scattering 
of  electrons,  which  has  already  been  shown8,9  to 
lead  to  asymmetries  in  the  thermopower. 


Fig.  4.  In  our  device,  the  well-established  magnetic- 
field  symmetries  under  constant  electrical  bias  are 
clearly  evident,  panel  (a).  The  newly  observed 
magnetic-field  symmetries  under  constant  thermal 
bias  seen  in  this  work  show  that  the  symmetries 
described  by  Safi  hold  true.  Panels  (b)  and  (d)  show 
expected  symmetries.  Panel  (c)  shows  two  traces  that 
are  not  expected  to  be  symmetric.  The  green  and  red 
curves  in  each  panel  are  offset  for  clarity.  Figure  taken 
from  Ref.  7. 
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D.5.  Magnetic-Field  Symmetries  at  Finite  Thermal  Bias 

To  further  study  the  thermoelectric  transport  symmetries  in  these  devices,  we  looked  at  the 
magnetic-field  symmetries  under  finite  thermal  biases,  i.e.  the  non-linear  regime.  To  quantify  the 
degree  of  symmetry  as  a  function  of  thermal  bias,  we  introduced  the  symmetry  parameter 

B 

where 

Note  that  Sap  ranges  from  -1  to  +1,  were  -1  represents  complete  anti-symmetry,  and  +1 
represents  complete  symmetry.  The  results  of  our  symmetry  analysis  are  shown  in  Fig.  5.  We 
observe  that  the  dominant  trend  of  Sa/3  is  to  decrease  with  increasing  thermal  bias.  However,  this 
trend  is  not  universal  to  all  of  the  symmetry  parameter’s  components;  a  few  of  the  off-diagonal 
components  remain  unchanged  or  even  show  a  slight  increase  in  symmetry  with  increasing 
thermal  bias.  The  origin  of  the  decreasing  symmetries  could  be  due  to  non-linear  thermoelectric 
effects,  which  we  estimate  to  begin  at  Vh  ~  1  mV;  however,  they  could  also  be  due  to  other 
effects,  such  as  heating  of  the  non-heated  terminals. 


Fig.  5.  Here  we  show  the  symmetry  parameter,  Sa^,  as  a  function  of  the  applied  heating  voltage.  The  diagonal 

components,  a  =  /?,  show  a  clear  decrease  as  the  thermal  bias  is  increased.  The  same  decrease  is  only  partially  seen  in 

the  off-diagonal  components,  a  ^  p. 
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E.  Conclusions. Our  results  have  shown  that  the  transverse  thermopower  of  the  four-terminal  junction 
with  a  central  scatterer  is  very  sensitive  to  the  symmetry  of  the  four  QPC  transmissions  functions.  We 
have  shown  that  as  a  first  order  conceptual  model,  the  four-terminal  junction  can  be  thought  of  as  two 
pairs  of  QPCs,  one  left  and  one  right.  This  model  is  capable  of  describing  both  the  observed  sign  change 
in  the  transverse  thermopower  upon  reversal  of  the  thermal  gradient,  as  well  as  the  observed  nonlinear 
effects  in  the  transverse  thermopowers.  To  fully  take  all  four  terminals  into  account,  we  have  put  forth 

a  multi-terminal  thermoelectric  scattering  theory,  which  we  have  used  to  extract  the  applied  temperature 
increase  as  well  as  the  effective  transverse  thermopower.  Finally,  we  have  experimentally  demonstrated 
the  fundamental  magnetic-field  symmetries  of  thermoelectric  transport  coefficients  under  thermal  and 
electrical  biases.  We  have  also  shown  that  these  symmetries  tend  to  decrease  with  increasing  thermal 
bias.  From  our  experiments  we  cannot  conclude  whether  this  symmetry  breaking  is  due  to  a  heating 
effect  (for  example  the  destruction  of  phase  coherence  due  to  inelastic  scattering)  or  due  to  a  nonlinear 
effect. 

F.  Ongoing  and  Future  Work 

F.l.  Second  Generation  Devices 

In  response  to  the  measurements  reported  in 
Section  D.l-3  and  Ref.  1,  we  fabricated  a  second 
generation  of  devices  designed  to  investigate  the 
symmetry  of  the  four-terminal  junctions,  Fig.  1(c), 
a  device  with  a  vertically  displaced  circular 
scatterer,  Fig.  6(b),  and  arrays  of  triangular 
scatterers.  The  symmetry  measurements  have 
been  completed  and  are  stated  in  Section  D.4-5 

The  circular  device  will  test  our  hypothesis  that 
the  exact  shape  of  the  scatterer  is  not  critical  to 
generating  the  observed  thermoelectric  effects. 

The  multi-terminal  scattering  theory  we  have 
developed,  and  the  electronic  and  thermal  bias 
measurements  on  device  2  support  this 
hypothesis.  In  particular,  for  the  said 
measurements,  we  have  seen  comparable 
transverse  thermovoltages  while  individually 
heating  each  of  the  four  terminals  in  the  device 
shown  in  Fig.  1(c).  We  also  have  preliminary 
electronic  measurements  on  the  circular  scatterer 
device  that  suggest  this  is  true,  Fig.  7. 

Devices  consisting  of  arrays  of  triangular 
scatterers  were  also  fabricated.  In  a  possible 
continuation  of  this  project,  these  devices  will  be 
used  to  explore  transverse  thermoelectric  effects 
in  a  scaled-up  functional  nanomaterial,  Fig.  6(b), 
as  well  as  whether  asymmetric,  nanoscale 
obstacles  can  be  used  to  partially  rectify  the  flow 
of  heat  carried  by  electrons,  Fig.  6(c). 
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Fig.  6.  Three  of  seven  new  device  structures  fabricated  to 
test  for  the  rectification  of  heat  flow,  (a)  Device  with  a 
circular  central  scatterer.  (b)  Array  of  triangular  scatterers 
used  to  test  functional  nanomaterial  properties,  (c)  Device 
structure  for  studying  the  rectification  of  heat  flow. 
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Applied  Lateral  Voltage  (mV) 

Fig.  7.  The  observed  transverse  voltage  response,  V24,  of  the  displaced  circular  scatterer  device,  Fig.  6(a),  due  to  a 
lateral  electronic  excitation  suggests  that  the  exact  geometrical  symmetry  of  the  scatterer  is  not  critical  to  generating 
the  observed  effects.  The  key  point  here  is  that  the  magnitude  of  the  observed  transverse  voltage  is  comparable  to 

that  from  the  device  with  a  triangular  scatterer. 


F.2.  Theoretical  Analysis 

We  are  collaborating  with  two  theorist  groups:  David  Sanchez,  Mallorca  and  Peter  Samuelsson, 
Lund.  Through  these  collaborations,  we  are  refining  the  multi-terminal  thermoelectric  scattering 
theory  described  in  section  D.3.  Using  this  model  and  the  measurements  discussed  in  Section 
D.3-4.,  we  are  actively  comparing  the  observed  experimental  symmetries  with  theoretical 
predictions. 

F. 3.  Numerical  Simulations 

Peter  Samuelsson’ s  group  is  conducting  numerical  simulations  of  the  symmetry  device  shown  in 
Fig.  1(c).  Their  simulations  use  a  recursive  Green’s  function  algorithm  to  model  the  scattering 
matrix  describing  the  junction.  Due  to  the  sensitivity  of  the  actual  device  to  its  dimensions  and 
stray  charges,  we  do  not  expect  to  achieve  complete  agreement  between  experiments  and 
simulations.  Instead,  we  intend  to  compare  the  statistical  properties  of  the  junction’s  transport 
properties,  such  as  the  magnitude  of  the  transverse  voltage  as  a  function  of  an  applied  magnetic 
field. 
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